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@ Precede et dispositif de recalage d'images utilisant une corrdlation de phase avec ponddration pour 
determiner un decalage. 
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@ Le proc^d^ de recalage selon I'invention met 
en oeuvre une correlation de phase, c'est-d-dire 
une correlation entre les transformees discretes 
^ deux dimensions d'une image de reference et 
d'une image courante normalisee par rapport a 
son module, puis une transformation de Fourier 
inverse de la matrice de conrelation obtenue 
pour trouver un pic de correlation dans le 
domaine spatial des images, dans lequel, pour 
diminuer les effets du bruit dans ce traitement 
et renforcer les composantes moyennes, une 
fonction de ponderation est appliquee a la 
matrice de correlation dans le plan transforme 
avant d'effectuer la transformation de Fourier 
inverse. 

L'invention s'applique au recalage en transla- 
tion d'images notamment d celles obtenues par 
des capteurs en infirarouges. 
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L'invention se rapporte au traitement d'images et a plus particulierement pour objet un precede et le dis- 
positif correspondant de recalage d'images utilisant la correlation de phase avec ponderation pour determiner 
un d^calage. 

Les precedes de recalage d'images connus sont de pres ou de loin bases sur une analyse globale des 

5 contours apparaissant dans les images et visent a determiner un decalage general entre contours pour remettre 
en coincidence ces contours associes aux memes objets de fagon a recaler une image courante par rapport 
a une image de reference et a detecter ainsi dans I'image recal^e les modifications eventuelles intervenues 
entre les deux images. Cette technique s'applique notamment aux traitements d'images obtenues d partir de 
capteurs infrarouges dans des applications a la surveillance parexemple. 

10 Les methodes les plus classiques de recalage d'images effectuent des traitements de correlation entre les 

deux images mais ce type de m^thode fait apparaTtre sur certaines des frequences des pics de correlation qui 
sont fonction de I'energie. On connait done une methode plus perfectionnee dite "methode par correlation de 
phase" qui est une methode de recalage utilisant des transformations de Fourier, un peu semblable aux metho- 
des de correlation classiques a la grosse difference pres que chaque composante est normalis6e par rapport 

15 d son module dans le domaine transforms de Fourier. Ainsi une composante de frequence quelconque a la 
meme importance que les autres pour Tobtention des pics de correlation, quelles que soient les amplitudes 
respectives de ces composantes de frequences. Un tel precede de recalage d'images dit par correlation de 
phase est decrit notamment dans le brevet frangais n"" 88 14675 au nom de la Demanderesse. 

Ce type de procSde de recalage a des performances convenables lorsque les images ne sont pas bruitSes. 

20 Par centre, ce precede perd de sen efFicacit6 en presence d'images de faible definition et fortement bruitSes 
(par exemple images infrarouge). 

L'invention a done pour objet un precede de recalage d'images mettant en oeuvre une correlation de phase, 
mais adapts pour rendre plus robuste la technique de recalage. Pourcela, l'invention effectue une pondSration 
du spectre de correlation de fagon ^ accroTtre influence de la zone spectrale susceptible de contenir les 

25 contours tout en reduisant I'lnfluence du bruit. Ainsi la methode de recalage cerrespondante est applicable aux 
images tres bruitees ou n'ayant que tres peu de definition, sans toutefois mettre en oeuvre des operations 
compliquees de prS-traitement avant la correlation. 

Selen l'invention, un precede de recalage d'images selen lequel une transformation de Fourier discrete a 
deux dimensions Fr(u,v) des valours de luminance fr(x,y) d'une image de reference ayant ete effectuee au prea- 

30 lable, le procede comporte : 

- une etape durant laquelle une transformee de Fourier discrete d deux dimensions Fc(u,v) des valeurs 
de luminance fc(x,y) de I'image courante est effectuSe 

- puis une etape de calcul de la matrice de correlation C(u,v) entre les transformees de Fourier de I'image 
de reference et de I'image courante, normalisees par rapport aux modules des coefficients obtenus 

35 - puis une etape de calcul de la transformee de Fourier inverse de la mafrice de correlation faisant appa- 

raTtre un pic de correlation dans le domaine spatial (x,y) caracteristique d'un deplacement entre elements 
d'images analogues, estcaractehse en ce que, avant I'etape de calcul de la transformee de Fourier inverse, 
la matrice de correlation est ponderee par une fonction de ponderation P(u,v) pour limiter Tinfiuence du 
bruit et renforcer les valeurs meyennes dans lesquelles se treuvent la plupart des contours. 

40 L'invention a egalement pour objet un dispositif de recalage d'images destine a la mise en oeuvre de ce 

procede. 

L'invention sera mieux comprise et d'autres caracteristiques apparaTtront a I'aide de la description qui suit 
en reference aux figures annexees. 

- La figure 1 est un schema synoptique du procede de recalage d'images selon ['invention ; 

45 - Les figures 2a, 2b et 2c sont des schemas explicatifs de I'etablissement de la fonction de ponderation ; 

- La figure 3 est un schema representant en trois dimensions la fonction de ponderation appliquee k la 
fonction de correlation dans le plan transforms ; 

- La figure 4 est le meme schema a une echelle differente et pour 1/4 seulement du plan transforms. 
Comme indiquS ci-dessus l'invention dScrite ci-apres met en oeuvre une version perfectionnSe de la tech- 

50 nique de corrSlation de phase a deux dimensions, particulierement destinSe aux images de mauvaise qualitS. 
Lorsqu'une mSthode de recalage opere de fagon globale sur les niveaux de gris, ce sont les contours meme 
de faible amplitude qui permettent d'aboutir au dSplacement vrai. Au moment ou s'effectue I'opSration de conrS- 
lation, il est intSressant d'une part de donner plus de poids aux zones de I'image qui sont supposSes contenir 
ces contours, et d'autre part de rSduire I'lnfluence du bruit. 

55 Dans la mesure ou Ton procede par corrSlation de phase, le spectre de la conrSlation est facilement acces- 

sible. Cette approche est done mise a profit dans la prSsente invention ou une pondSration du spectre de la 
corrSlation est effectuSe juste avant sa transfomnation de Fourier inverse. La description dStaillSe qui suit va 
permettre de mieux expliquer l'invention. 
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Si Ton se refere a la figure 1 qui illustre le schema synoptique de la correlation avec la ponderation du spec- 
tre de correlation selon {'invention, il est possible dans un premier temps de reperer les etapes qui existaient 
d^jd selon la m^thode classique de recalage utilisant une conr^latlon de phase pour determiner un d^caiage 
en translation apres, avoir au prealable determine un angle de rotation. On pourra se referer a la demande de 

5 brevet frangais precitee n° 88 14675 pour une description detaillee de cette methode incorporee a la presente 
description par reference. Ces etapes sont rappelees ci-apres. 

A partir d'une memoire d'image 10 dans laquelle est memorisee Timage dite courante, c'est-a-dire Timage 
d recaler par rapport S une image de reference, le proc^d^ consiste dans une premiere ^tape 11 ^ appliquer 
la transfomnee de Fourier discrete a deux dimensions a I'image de reference d'une part et a I'image a recaler 

10 d'autre part. Ces images etant supposees constituees de NxM points, chaque coefficient transforme note 
Fr(u,v) pour I'image de reference est ^gal d : 

M-1 N-1 

V 



15 



Fr(u,v) = l/N.l/M 2_ fr{x,y).Exp(-21Ij,(xu/N * y v/M) ) 

x=0 y=0 



avec u = 0 N-1 

etv = 0 M-1 

20 et ou fr(x,y) est la luminance du pixel de coordonnees (x,y) de I'image de reference. 

De meme pour I'image courante, la transformee de Fourier de cette image courante notee Fc(u,v) est egale 

25 M-1 N-1 

Fc(u,v) = l/NA/uY2 Yl fc(x,y).Exp(-21Ij.(xu/N + yv/M)) 
x=0 y=0 

30 OU fc(x,y) est la luminance du pixel de coordonnees (x,y) de I'image courante. 

A partir de ces deux transformees de Fourier, une 6tape 12 consiste k calculer rintercorr^lation entre ces 
deux fonctlons de la mani^re suivante : 

C(u,v) = (Fr(u,v).Fc*(u,v))/|Fr(u,v).Fc(u,v)|. 
ou : Fc*(u,v) est le complexe conjugue de Fc(u,v). 
35 Comme indiqu^ cl-dessus, cette fonctlon de correlation dite "de phase" est normalis^e au niveau de ses 

composantes par rapport aux modules de ces dernieres. 

L'etape suivante 13 consiste alors a calculer la transformee de Fourier inverse a deux dimensions de la 
fonction d'intercorrelation C(u,v) telle que : 

c(x,y) = F-i(C(u,v)) 

40 pour revenir au domaine spatial (x,y), c(x,y) etant le plan de correlation contenant le ou les pics caracteristiques 
des decalages entre les deux images. La derniere etape 14 du precede classique de correlation consiste a 
estimer un vecteur de decalage de composantes dx et dy entre I'image de reference et I'image courante dans 
le plan de correlation c(x,y). II est possible de noter que rapplication d'une fenStre de Hamming sur les deux 
images, image de reference et image courante, avant de calculer les transformations de Fourier discretes peut 

45 avantageusement etre effectuee. 

Les etapes decrites ci-dessus sont reproduites dans le procede selon I'invention mais k partir de la fonction 
d'intercorrelation C(u,v) des transformees de Fourier discretes a deux dimensions de I'image de reference, 
Fr(u,v), et de I'image courante, Fc(u,v), une ponderation de I'intercorrelation est effectuee dans une etape 15 
par une fonction de ponderation P(u,v) a deux dimensions appliquee a I'intercorrelation C(u,v) de la maniere 

50 suivante : 

C'(u,v) = C(u,v).P(u,v) 

L'interconrelation dans le domaine spatial est alors calculee par la transformee de Fourier inverse a deux 
dimensions de I'intercorrelation ponderee dans le domaine de Fourier C'(u,v), le vecteur de decalage etant 
extrait de Tintercorreiation dans le domaine spatial de la meme maniere que dans Tart connu. 
55 La fonction de ponderation appliquee a I'intercorrelation est optimisee pour obtenir simultanement une 

reduction de influence du bruit et un renforcement des composantes a frequences moyennes qui sont celles 
dans lesquelles apparaissent la plupart des contours. Pour cela, d titre d'exemple la ponderation peut etre cal- 
culee de la maniere suivante : 
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P(u,v) est construite a partir de fonctions de transfert de type filtres de Butterworth caracterisant res- 
pectivement un filtre passe-bas et un filtre passe-haut, de pseudo-frequences de coupure convenables, une 
valeur r^siduelle pour les basses frequences 6tant conserv^e sous la forme d'un coefficient constant ajout^ k 
la caracteristique du filtre passe-haut. 
5 Les courbes representees sur les figures 2a, 2b et 2c representent respectivement, pour une dinnension, 

les formes des courbes de reponse des filtres avec les pseudo-frequences de coupure representees par les 
coefficients Di et D2 et la fonction de ponderation resultant du produit de ces deux caracteristiques toujours 
sur les mSmes bases. 

Sur la figure 2a a ete represents la courbe Pi(u) du filtre passe-haut avec une pseudo-fr6quence de cou- 
10 pure egale a Di et exprimee dans le plan (u,v) a deux dimensions par 

Pi(u,v) = 0,25 + 0,75W(1 +0,414.Di/D(u,v))2N1 

ou D(u,v) = Vu2 + v2 et ou Ni est un nombre. 
A titre d'exemple Di = 5 et Ni= 2. 
^5 Sur la figure 2b a ete representee la courbe P2(u) du filtre passe-bas avec une pseudo-frequence de cou- 

pure 6gale d D2 et exprimee dans le plan (u,v) ^ deux dimensions par : 

P2(u,v) = 1 W( + 0,414.D(u,v)/D2)2N2 

ou N2 est un nombre. 

A titre d'exemple D2 = 35 et N2 = 5. 

20 La figure 3 illustre les coefficients de ponderation dans le plan transfonme pour une transformee de taille 

128x128 pointsja figure 4 illustrant la meme fonction de ponderation suivant 1/4 du plan transfonme. Dans 
I'exemple illustr^ la zone spectrale permettant de privil^gier la correlation est comprise entre Di = 5 et D2 = 
35 ; mais d'autres valeurs peuvent etre choisies dependant du type de capteurs d'ou sont issues les images, 
voire du contenu de I'image elle-meme. Cette precision est donnee car il faut veiller lors du traitement a ne 

25 pas supprimer des zones spectrales qui sont susceptibles de contenir rinfomriation n^cessaire ^ I'estimation 
du ou des vecteurs de decalage. 

II resulte des explications indiquees ci-dessus que la fonction de ponderation se presente pour le traitement 
sous la forme d'un tableau de nombres reels utilises comme coefficients de la fonction de correlation appliques 
aux differents points du plan de correlation en deux dimensions, ce tableau de nombres reels etant tel qu'en 

30 presence de transform^es de tallies NxM, la taille de la m^moire de stockage des coefficients pourra Stre de : 

- N.(M/2+1) en tenant compte de la symetrie hermitienne ; 

- (N/2+1). (M/2+1) si la fonction d'intercorrelation est de plus de symetrie circulaire ; 

- (N/2+1).(M/4+1) avec une fonction de ponderation de symetrie circulaire. 

II apparaTt clairement de la description ci-dessus que la fonction de ponderation appliquee a la fonction de 
35 rinter-correiation k deux dimensions decrite ci-dessus n'est qu'un exemple de realisation et peut etre adaptee, 
au contenu des images a traitor notamment. 

Le dispositif de mise en oeuvre du precede de recalage d'image selon Tinvention est analogue a celui qui 
a ete decrit dans la demande de brevet firangais n"" 88 14675 au nom de la Demanderesse a ceci pres qu'il 
comporte en outre une memoire de stockage des coefficients de ponderation telle que decrite ci-dessus. 

40 

Revendications 

1. Precede de recalage d'images selon lequel une transformation de Fourier discrete a deux dimensions 
^5 Fr(u,v) des valeurs de luminance fr(x,y) d'une image de reference ayant ete effectuee au prealable, le pro- 

cede comporte : 

- une etape durant laquelle une transfomriee de Fourier discrete d deux dimensions Fc(u,v) des valeurs 
de luminance fc(x,y) de I'image courante est effectuee 

- puis une etape de calcul de la matrice de correlation C(u,v) entre les transformees de Fourier de 
so I'image de reference et de Timage courante normalisee par rapport aux modules des coefficients obte- 

nus 

- puis une etape de calcul de la transformee de Fourier inverse de la matrice de correlation faisantappa- 
raltre un pic de correlation dans le domaine spatial (x,y) caracteristique d'un deplacement entre ele- 
ments d'images analogues, 

55 caracterise en ce que, avant I'etape de calcul de la transformee de Fourier inverse, la matrice 

de correlation est ponderee par une fonction de ponderation P(u,v) pour limiter I'influence du bruit et 
renforcer les valeurs moyennes dans lesquelles se trouvent I'information de structure. 
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Precede selon la revendication 1 , caracterise en ce que la fonction de ponderation dans le plan transforme 
P(u,v) est telle qu'elle reduit I'influence des connposantes de la nnatrice de correlation associees a des 
nombres de points faibles, inf^rieurs k une premiere valeur Di d'une mani^re analogue k un filtrage passe- 
haut, et reduit tres fortement I'influence des composantes de la matrice de correlation associees a des 
nombres de points eleves, superieur a une seconde valeur D2, d'une nnaniere analogue a une filtrage 
passe-bas, les composantes de la matrice associees a des nombres de points compris entre les deux 
parametres Di et D2 etant privilegiees. 

Precede selon la revendication 2, caracterise en ce que les parametres de la fonction de pond6ration sont 
determines a partir des caracteristiques du capteur dont sont issues les images a traitor. 

Dispositif de recalage d'images destine a la mise en oeuvre du precede selon I'une des revendications 1 
a 3, caracterise en ce qu'il comporte une memoire PROM de stockage des coefTicients de ponderation 
formant la fonction de ponderation P(u,v). 

Dispositif selon la revendication 4, caracterise en ce que, pour des images comportant N.M points, la taille 
de la matrice est limit^e ^ N.(M/2+1), en tenant compte de la sym^trie hermitienne. 

Dispositif selon la revendication 5, caracterise en ce que la taille de la matrice est limitee k (N/2+1 ). (M/2+1) 
en tenant compte de la symetrie hermitienne et lorsque la fonction de ponderation est de symetrie circu- 
laire. 

Dispositif selon la revendication 6, caracterise en ce que la taille de la matrice est limitee a (N/2+1 ). (M/4+1) 
en tenant compte de la symmetric hermitienne, lorsque la fonction de ponderation est de symetrie circu- 
laire. 
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